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Abstract 

Non-commutative black holes are characterised by a minimum mass which would result in a 
remnant after the Hawking evaporation ends. We numerically study the decay of neutral non- 
commutative black holes for up to ten spatial dimensions and typical parameters that would make 
. their production possible at the LHC. Neglecting possible accretion mechanism, we find that decay- 

times are extremely short. 

o . 

^ ■ 1 Introduction 

One of the most striking predictions of General Relativity is the existence of black holes, that is event 
horizons which surround and hide the point-like or otherwise singular classical sources. On general 
quantum-mechanical arguments, one however expects that singular sources are smeared out and there 
are indeed many theoretical reasons to believe the fabric of space-time can be effectively described by 
non-commutative geometry at short scales (see, for example Section 2 of Ref. pQ and Ref. [2]). Such 
modifications to General Relativity must necessarily have an impact on the structure of black holes. 

The non-commutative analogue of the Schwarzschild black hole metric was found in Ref. [3] (see also 
Ref. [3]) and subsequently generalised to include the electric charge [5] and extra-spatial dimensions [51 [7]. 
A common feature for all such metrics is that there exists a minimum mass for the black hole which 
solely depends on the non-commutative length ^0 and the fundamental scale of gravity M g — £~ 1 (for a 
comprehensive and updated review, see Ref. [1]). The former is usually identified with the Planck length 
£ p and the latter with the Planck mass M p in four space-time dimensions, but £ could actually be much 
larger and M g much lower, possibly around 1 TeV, if extra-dimensions exists [5J [S]. This has opened up 
the possibility of having micro-black holes (mBH) with a mass of a few TeV's QIJj that might be produced 
and detected at the Large Hadron Collider (LHC) [TTl [T2l [T5] . 

The standard picture of mBH at the LHC is that they should be produced with a mass of around 
10 TeV and very rapidly decay completely [11] via the Hawking effect [14]. Generally small (but possibly 
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significant) corrections to this scenario arise from the use of the microcanonical picture [15] , which 
essentially enforces energy conservation throughout the evaporation process and suppresses the emission 
for small mBH mass |12[ 116] . If the non-commutative description of mBH given in Ref. pQ is correct, 
one however expects that a (stable) remnant will be left at the end of the Hawking evaporation [T?] . 
which might exit the accelerator and detectors. In this context, and for phenomenological purposes, it is 
therefore important to estimate the typical decay-times for such objects. For example, the issue of mBH 
life-times has recently been debated in Refs. [18] . 

Since it is widely accepted that mBH discharge very quickly, via several processes, including the 
Schwinger mechanism at least for 3 < d < 5 [7], we shall here consider only the neutral metrics of 
Refs. [U H] for d = 5, ... , 10 spatial dimensions (noting that higher values of d seem to be favoured [19]). 
In order to have mBH produced at the LHC, we need a minimum mBH mass of order 1 TeV, which 
in turn implies that l g ~ I. Further, we shall see that the typical temperature of such mBH's always 
remains smaller than their mass, and microcanonical corrections can thus be neglected. The evolution of 
the mBH mass is then determined by solving the relevant equation numerically for initial mBH mass of 
order 10 TeV and an upper bound for the decay-time obtained. 

We shall use units with c = h = 1 . 



2 Non- commutative neutral mBH 

We recall here that the non-commutative metric of a neutral black hole of (asymptotic) proper mass m 
in d + 1 space-time dimensions is given by [7J 

ds 2 = -A d dt 2 + A' 1 dr 2 + r*- 1 d^_ x , (2.1) 

where 

2G d m (d r 2 \ . . 



r d - 2 r(d/2) 1 \r i 2 

and 



7(0; a;) = / u^e^du . (2.3) 



With d = 3 spatial dimensions, the Newton constant G3 = Gn = £ P /M p = £ 2 and one usually assumes 
that £ ~ £ p . In the presence of extra spatial dimensions, d > 3, we can likewise write Gd = £ g ~ 2 '/M g = 
£g X , but keeping £ and £ g distinct. Of course, the non-commutative length must be short enough to 
agree with present experimental bounds on the validity of the Newton law at short distance, that is 
£ < 1 /im [20] and for the gravitational mass we shall assume M g ~ 1 TeV. 
The radial mass-energy function for the above solution is given by 

, . m / d r 2 \ , 

m ^=wmA*p) ' (2 - 4) 

and 

lim m(r) — m , (2-5) 

r — >oo 

as anticipated. Further, for sufficiently large proper mass, there exists two horizons, the outer one being 
located at r = r + , which is related to m by 



Ma - £ d - 2 - - r d - 2 F (d/2) ('>(,) 
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where T(o) denotes the Gamma function and Md = Gdm the Arnowitt-Deser-Misner (ADM) mass of the 
black hole 0. The inverse Hawking temperature (in geometrical units) is finally given by 



4 7T T a 



(*-2)-2(^y 



(2.7) 



3 Hawking evaporation 

In order to study the evaporation of a mBH, we need to estimate the corresponding Hawking flux that 
is the amount of energy emitted through a unit area per unit time in the form of thermal radiation. The 
latter is described by quantum fields propagating on the (d + l)-dimensional background (|2.ip and we 
recall that the space-time non-commutativity makes the field theory UV finite, thanks to the presence 
of a damping term in the momentum-space propagator, namely G(p) ~ exp f— £ 2 p 2 /8) (for the details, 
see Ref. [H]). As a result, considering for simplicity the case of a massless scalar field and neglecting the 
grey-body factor, one finds 



3> = 2 / 1 % e ~" e2p2p ■ (3.1) 
J (2 7r) d v ' 

From this expression, one can then obtain the total luminosity Cd which governs the time dependence of 
the mBH proper mass in the canonical picture [7], 

- — oc F d (r + ) [j^ J r(2n + 4)C(2n + 4)^/: d (r + ), (3.2) 

where £(a) denotes the Riemann zeta function. 

The exact form of the function Fd in the above expression is actually difficult to determine and for 
a number of reasons. To begin with, one should consider the contribution of all the different kinds of 
particles that can be emitted, like Standard Model fermions and bosons as well as (bulk and brane) gravi- 
tons. Further, particles with different spins propagate differently and their emission is thus suppressed by 
spin-dependent grey-body factors, which, in turn, strongly depend on the precise near-horizon geometry. 
In this respect, there are arguments which suggest that Standard Model particles propagating only in our 
(3 + I )-dimensional space-time makes for most of the Hawking radiation [55], so that F d oc (the area 
of the brane-section of the horizon), as well as arguments against this scenario [53], according to which 
one should instead have Fd oc (the area of the bulk horizon, neglecting a possible squeezing [24]). 
Finally, one realistically expects that mBH at the LHC would be produced with an intrinsic spin (see, 
for example, Ref. [25]), but no rotating analogue of the metric (|2.1|) is known to date. 

Since our main goal here is to provide order of magnitude estimates of the typical decay-times, we 
shall consider three possibilities: 

B) The mBH emits all the particles in the entire bulk: 

Ff=r^(3-^. (3.3) 

b) The mBH emits like a four-dimensional black hole except for the modified (inverse) temperature: 

F h d =rl(3-/. (3.4) 

bB) The mBH evaporates with a dependence on the temperature like that in Eq. ()3.3|1 and on the 
2-dimensional horizon area like in Eq. (|3.4p : 



2 Note that Md has dimensions of (length)" 
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All the above forms will be used for the numerical simulations reported on later. 

In order to solve (|3.2p numerically, we introduce dimensionless quantities (denoted by a bar) by 
expressing everything in units of i (to a suitable power). For example, the dimensionless ADM mass will 
be given by 

M d = M d £ 2 - d , (3.6) 

and the proper mass by 

m = m£ = M d (J-^ . (3.7) 

We shall sometimes find it more convenient to express the proper mass m in units of M g (denoted by a 
tilde) as 

*-£-*(£)"• (3 - 8) 

Since we assume M g ~ 1 TeV, rh will just be the mBH mass in TeV's. Next, we shall need to truncate the 
series in Eq. (|3.2[) to a maximum integer value n = n max . With the above assumptions and redefinitions, 
Eq. (|3~2|) becomes 

Ajf, "max / -, \ n 

= -cF d (f + ) ]T V^r^" 2nr(2n + 4)c(2?1 + 4) = - ld{f+) ' (3 - 9) 

ra=0 



where C is a numerical constant which we shall comment upon later. Finally, by making use of Eq. (|2.6p . 
we can obtain an equation which contains the horizon radius as the only time-dependent variable, 



df+ _ df+ Pdi{i,f 2 +) ft 



d-l 



at dm 2ixr\ l T(d/2) \ t ) 

which we now proceed to study numerically, with the aim of determining the evolution of a mBH with 
initial mass m(t = 0) ~ 10 TeV. 



4 Numerical results 

Upon employing the three forms listed in the previous Section for the function F d , we always found 
the same qualitative picture. In particular, for fixed C, the choice of in Eq. (|3 . 3|) always produced 
the longest decay-times for 5 < d < 10. Moreover, for a given initial mBH mass and choice of F d , the 
decay-time increases for increasing d and is (roughly) linearly proportional to C. Since larger values of 
d are anyway favoured [19j, we shall therefore show the detailed analysis just for the case d = 10 and 
Fiq = F® Q with C = 1. This choice yields the longest decay-time overall (for fixed C = 1), which can 
thus serve as a "worst" (or "best") case scenario, and we shall further comment about other cases when 
relevant. 

First of all, we need the minimum mBH mass. This is obtained by first minimising the ADM mass (|2.6p 
with respect to r+, which, for d = 10, yields the minimum (dimensionless) horizon radius f m j n ~ 1.11, 
corresponding to a minimum ADM mass M m in — 133. In order for the latter to translate into a minimum 
proper mass m m in — 1.2 TeV, we thus set (see Eq. (j3.8JI ) 

e g ~1.8£. (4.1) 

For smaller values of d, one obtains slightly smaller ig's, but with a very weak dependence on the space 
dimension. In fact, for d = 5, one has the smallest value £ g ~ 1AI. 
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Figure 1: Temperature (left panel) and its ratio with proper mass (right panel) vs horizon radius in 
d = 10 (where f + > 1.11). The maximum temperature ~ 0.11m is reached for f + ~ 1.68. 



In Fig. [TJ we then plot the mBH temperature /J^ 1 (left panel) and its ratio with the proper mass fh 
as functions of the horizon radius f + (right panel). Since the latter ratio becomes smaller for lower d, it 
is clear that the mBH mass always remains significantly larger than the typical energy of emitted quanta 
(proportional to the temperature) and the canonical expression (|3.2[) applies throughout for all values of 
d (and all choices of F d ). 

Next, we need to determine how many terms in the sum in Eq. (|3.10[) we should keep. In Fig. [21 left 
panel, we show the dimensionless luminosity in d = 10 with Fiq — F^ and C = 1 for n max = 0,2 and 10 
around the peak, where discrepancies are the largest. We can thus conclude that it is well sufficient to 
use n max = 2 (the relative difference with respect to n max = 10 does not depend on the choice of Fd and 
C and is less than 8 • 10 -3 ; see right panel). For smaller values of d, this approximation actually becomes 
even better. For example, for d = 5 one obtains the smallest relative error of less than 6 • 10~ 5 . 

From Eq. (|2.6p and (|3.8[) . we find that an initial mBH proper mass of about lOTeV corresponds to 
an horizon radius f + ~ 2.5 for d = 10. With this initial condition, Eq. (|3.10|) with C = 1 and Fio = F^ Q 
yields the solution plotted in the left panel of Fig. [3l from which the time-dependent proper mass can 
be easily obtained (right panel). The latter plot shows that the decay occurs through several stages: it 
is initially rather slow (for < t < 300), subsequently becomes much faster (for 300 < t < 10 6 ) and 
finally approaches the remnant again slowly (for t > 10 6 ). For our choice of £ g — ITeV -1 ~ 1.8 
we obtain that i ~ 0.6 TeV~ x ~ 3 • 10 _28 sec and the estimated decay-time ~ 10 10 ^ ~ 10 _18 sec is 
practically instantaneous. This result changes rather weakly with d, decreasing down to the minimum 
<5 ~ 10 6 £ ~ 10~ 22 sec. As we mentioned before, different choices of Fd produce slightly different decay 
curves, but do not yield longer decay-times. For example, for d = 5 and F5 = F§ in Eq. (|3.4|) one obtains 
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Figure 3: Left panel: horizon radius for initial mBH mass m ~ 9.2 TeV (solid line) and its minimum 
value (dotted line). Right panel: mBH proper mass for same initial condition (solid line) and remnant 
mass m min ~ 1.2 TeV (dotted line). 

the shortest possible value for C = 1, that is t$ ~ 10 4 1 ~ 10~ 24 sec. For different values of C, all the 
decay-times scale roughly linearly. Assuming C > 10~ 2 and considering the case d = 10 with Fio = F® Q 
thus yields the fairly conservative upper bound 

idocay < KT 16 sec , (4.2) 

which is extremely short compared to the sensitivity of present detectors (on the order of hundreds of 
picoseconds). 



5 Conclusions 

We have estimated the decay-times of non-commutative mBH's that might be produced at the LHC and 
found that they should evaporate nearly instantaneously (see the upper bound in Eq. I|4.2|) ). much the 
same as is expected according to the "canonical" scenario of Ref. [11]. We can consequently conclude 
that non-commutative mBHs would evaporate within the detectors, since they can propagate at most a 
few nanometers away from their point of production during the evaporation. 

The present results are similar to previous estimates of decay-times for usual ("commutative") mBH's 
in the ADD brane- world [8] , described according to the microcanonical picture [15] , for which life-times 
were obtained on the order of 10 _17 sec or shorter [12]. Remarkably, our conclusions are also compara- 
ble with the recently estimated life-times of mBH's derived from Generalised Uncertainty Principles in 
Ref. [26j . Actually, this is not surprising, since the dependence of the temperature on the mBH mass 
is roughly similar in all of these cases [IB]- Of course, the mBH's of Ref. [T2] evaporate completely 
(thus reaching zero temperature for vanishing mass) and their life-time is actually dominated by the 
latest stages (when m < 5 TeV), whereas the end-point of the process here, as well as in Ref. [26], is a 
(presumably) stable remnant of mass m min ~ 1 TeV and zero Hawking temperature. 

Let us conclude by mentioning that no accretion mechanism has been included in our analysis of the 
evaporation, which is fully consistent given the very short decay-times. What happens afterwards and 
the final fate of the remnant are still open questions. 
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